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1 Governing equations for waves on the sea surface 

In this chapter we shall model the water as an inviscid and incompressible fluid, and 
consider waves of infinitesimal amplitude so that the linearized approximation suffices. 

Recall in the first chapter that when compressibility is included the velocity potential 
defined by u = V$ is governed by the wave equation: 



where c = a/ 'dp/ dp is the speed of sound. Consider the ratio 

c 2 dt 2 CO /K 



V 2 $ c 2 

As will be shown later, the phase speed of the fastest wave is uj/k = sfgh where g is the 
gravitational acceleration and h the sea depth. Now h is at most 4000 m in the ocean, 
and the sound speed in water is c = 1400 m/sec 2 , so that the ratio above is at most 
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We therefore approximate (1.1) by 

V 2 $ = (1.2) 

Let the free surface be z = ((x,y,t). Then for a gently sloping free surface the 
vertical velocity of the fluid on the free surface must be equal to the vertical velocity of 
the surface itself, i.e., 

m = dz-> z = - (L3) 

Having to do with the velocity only, this is called the kinematic boundary condition. 
For small amplitude motion, the linearized momentum equation reads 

du 

p-^ = -VP-pge z (1.4) 
Now let the total pressure be split into static and dynamic parts 

P = p +P (1.5) 



where p is the static pressure 



which satisfies 



p = -pgz (1.6) 



= -Vp + -pge z (1.7) 

Let us assume that wind is not present, so that the air above the sea surface is essentially 
stagnant. The static pressure is hydrostatic, which can be taken to be zero at z = 
without loss of generality. It follows that 

du dV$ 

p m= p ^r = - Vp (L8) 

so that 

" = -"m < L9 > 

Because of its very small density the passive air motion due to surface waves has negible 
effect on wate motion. We assume the dynamic air pressure to be zero on the free 
surface. 

It will be shown shortly that for sufficiently long waves, surface tension can be 
ignored. Continuity of pressure requires that 

p = Q, z = (. 



to the leading order of approximation, we have, therefore 



P9( + p-j£ =0, z = 0. 



(1.10) 



Being a statement on forces, this is called the dynamic boundary condition. The two 
conditions (1.3) and (1-10) can be combined to give 



d 2 $ 5$ 



If surface tension is also included then we adopt the model where there is a thin 
film covering the water surface with tension T per unit length. Consider a horizontal 
rectangle dxdy on the free surface. The net vertical force on water from four sides of 
the small rectangle is 



T 



dx 



x+dx 



dx 



dy+[T 



dy 



T 



y+dy 



dx 



dx 



d 2 ( d 2 ( 
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dx 2 dy 2 



dx dy 



Continuity of vertical force on an unit area of the surface requires 

p + t(— + — 1 =o. 

\ dx 2 dy 2 ) 

Hence 



p9C ~ p m +T \dx 2 



dy< 



0, z = 0. 



which can be combined with the kinematic condition (1.3) to give 



d 2 $ 



9$ 



g- 



d 2 



d 2 \ a$ 



o. z = o 



(1.12) 



(1.13) 



dt 2 dz p \dx 2 dy 2 J dz 
When viscosity is neglected, the normal fluid velocity vanishes on the rigid seabed, 



n • V$ = 

Let the sea bed be z = —h(x, y) then the unit normal is 

_ (h x ,.h y ,l) 



(1.14) 



Hence 



^TThlTh 2 
d$ dh 9$ dh 9$ 



dz dx dx dy dy 



z = -h{x,y) 



(1.15) 



(1.16) 
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2 Progressive waves on a sea of constant depth 



2.1 The velocity potential 

Consider the simplest case of constant depth so that 

— = 0, z = -h. (2.1) 

and sinusoidal waves with infinitively long crests parallel to the y axis. The motion is 
in the vertical plane (x,z). Let us seek a solution representing a wavetrain advancing 
along the x direction with frequency uj and wave number k, 

$ = f( z )e lkx - luJt (2.2) 

In order to satisfy (1.2), (1.13) and (2.1) we need 

f" + k 2 f = 0, -h<z<0 (2.3) 

-co 2 f + gf' + -k 2 f' = 0, z = Q, (2.4) 
P 

/' = 0, z=-h (2.5) 
Clearly the solution to (2.3) and (2.5) is 

f(z) = B cosh k(z + h) 

implying 

$ = ficoBhjfc(2 + /i)e ite - <wi (2.6) 
In order to satisfy (2.4) we require 

co 2 = \gk + -k 3 ) tanh kh (2.7) 
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which is the dispersion relation between w and k. From (1.3) we get 



dt dz 



= {Bksmhkh)e lkx - iuJt {2.1 

z=0 



Upon integration, 

£ _ j^ e ikx-iu)t _ Bk sinh kh ^ ikx _ iu)t 
— «(J 
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where A denotes the surface wave amplitude, it follows that 

B ~ IUA 



k sinh kh 

and 

$ = ~ %u coshk(z + h)e ikx - ialt 
k smh kh 

-igA ( Tk 2 \ coshk(z + h) 



, 1 + ; ' > e lkx ~ luJt (2.10) 

u) \ gp J cosh kh 

2.2 The dispersion relation 

Let us first examine the dispersion relation (2.7), where three lengths are present : the 
depth h, the wavelength A = 2ir/k, and the length X m = 27r/k m with 

k m = M, \ m = ^L = 2irJ- (2.11) 
V T k m y gp 

For reference we note that on the air-water interface, T/p = 74cm 3 /,s 2 , g = 980 cm/s 2 , 
so that X m = 1.73cm. The depth of oceanographic interest ranges from O(10cm) to 
thousand of meters. The wavelength ranges from a few centimeters to hundreds of 
meters. 

Let us introduce 



9P 

u m = £gn m = *g\i 

then (2.7) is normalized to 



co 2 m = 2gk m = 2g x /^ (2.12) 



to 2 1 k ( k 2 \ 

— = -— 1 + — tanh^ (2.13) 

Consider first waves of length of the order of A m . For depths of oceanographic 
interest, h ^> A, or kh ^> 1, t&nhkh ~ 1. Hence 

or, in dimensional form. 



u m 2 k m \ k^ 



TP 

UJ 2 = gk+ — (2.15) 

p 



The phase velocity is 



u q ( Tk 2 . 
c=T = l 7 U + — ( 2 - 16 ) 



A; V k \ gp 
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Defining 

the preceding equation takes the normalized form 



Vn = ^ (2-17) 



Cm y 2 ( A: 

Clearly 

iTk 

c^W— , if k/k m y>>l, or A/A m <l (2.19) 

Thus for wavelengths much shorter than 1.7 cm, capillarity alone is important, These 
are called the capillary waves. On the other hand 

|, if k/k m < 1, or A/A m < 1 (2.20) 

Thus for wavelength much longer than 1.73 cm, gravity alone is important; these are 
called the gravity waves. Since in both limits, c becomes large, there must be a minimum 
for some intermediate k. From 

dc 2 g T 



the minimum c occurs when 



= = k m, or A = A m (2.21) 



The smallest value of c is c m . For the intermediate range where both capillarity and 
gravity are of comparable importance; the dispersion relation is plotted in figure (1). 
Next we consider longer gravity waves where the depth effects are essential. 



CO 



^gktanhkh (2.22) 



For gravity waves on deep water, kh 3> 1, tanh kh — > 1. Hence 

lo w ^/gk, c ~ ( 2 - 23 ) 

Thus longer waves travel faster. These are also called short gravity waves. If however 
the waves are very long or the depth very small so that kh 1, then tanh kh ~ kh and 



bj^akygh, ctzygh (2.24) 




Figure 1: Phase speed of capillary-gravity waves in water much deeper than A TO . 

Form intermediate values of kh, the phase speed decreases monotonically with increasing 
kh. All long waves with kh <C 1 travel at the same maximum speed limited by the depth, 
■\/gh, hence there are non-dispersive. The dispersion relation is plotted in figure (??). 



2.3 The flow field 

For arbitrary k/k m and kh, the velocities and dynamic pressure are easily found from 
the potential (2.10) as follows 

d$ _ gkA ( A _ Tk 2 \ coshk(z + h) Jkx _ iaJt 



u = — = ± 1 + ^_Lv e .te-wi ( 2 25) 

ox u \ gp J cosh kh 

d$ _ -igkA t Tk 2 \ smhk(z + h), kx _ lidt 



w = — = —2 1 + ^e !M - lut (2.26) 

\ gp J cosh fen 



v = - p m =P9A { 1 + ^ cosh kh e (2 - 27) 

Note that all these quantities decay monotonically in depth. 
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Figure 2: Phase speed of capillary-gravity waves in water of constant depth 

In deep water, kh ^> 1, 

u = ^ ( 1 + —) e kz e ikx -^ (2.28) 
u \ gp J 

w = ^ = zMM ^ + 1^ W--< (2 . 29) 

dz u \ gp J 

(9<3) / Tk 2 \ 
p = -p— = pgA 1 + e kz e ik X -ia>t (2.30) 

dt \ gp J 

All dynamical quantities diminish exponentially to zero as kz — )■ — oo. Thus the fluid 
motion is limited to the surface layer of depth 0(A). Gravity and capillary-gravity waves 
are therefore surface waves. 

For pure gravity waves in shallow water, T = and kh 1, we get 

9^A ikx _ iuJ t 



u = e lkx ~ luJt (2.31) 



w = 0, (2.32) 

P = -P^r = P9 A ^ X ~ ZWt = P9C (2-33) 
ot 

Note that the horizontal velocity is uniform in depth while the vertical velocity is neg- 
ligible. Thus the fluid motion is essentially horizontal. The total pressure 

P = Po +P = pg((-z) (2.34) 

is hydrostatic and increases linearly with depth from the free surface. 
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2.4 The particle orbit 

In fluid mechanics there are two ways of describing fluid motion. In the Lagrangian 
scheme, one follows the trajectory x, z of all fluid particles as functions of time. Each 
fluid particle is identified by its static or initial position x ,z . Therefore the instan- 
taneous position at time t depends parametrically on x ,z . In the Eulerian scheme, 
the fluid motion at any instant t is described by the velocity field at all fixed positions 
x, z. As the fluid moves, the point x, z is occupied by different fluid particles at different 
times. At a particular time t, a fluid particle originally at (x ,z ) arrives at x,z, hence 
its particle velocity must coincide with the fluid velocity there, 

— = u(x,z,t), — = w(x,zA) (2.35) 
dt v ' dt K 5 ' ' K ' 

Once u, w are known for all x, z, t, we can in principle integrate the above equations to 
get the particle trajectory. This Euler-Lagrange problem is in general very difficult. 

In small amplitude waves, the fluid particle oscillates about its mean or initial posi- 
tion by a small distance. Integration of (2.35) is relatively easy. Let 

x(x ,z ,t) = x + x'(x ,z ,t), andz (x o, z ,t) = z + x'(x ,z ,t) (2.36) 

then x' <^C x, z' <C z in general. Equation (2.35) can be approximated by 

= u(x , z , t), ^77 = w(x , z , t) (2.37) 

From (2.25) and (2.26), we get by integration, 

9kA / + TP\ coshk(z + h) cikXo _ iu>t 

—iu) 2 \ gp ) coshfc/i 

gkA ( Tk 2 \ cosh k(z + h) . 

' 1 H ) j—j-j. sm(kx — cot) (2.38) 

(2.39) 

^) r ikx —iu)t 

^ \ 9P J cosh kh 

qkA ( Tk 2 \ sinh k(z + h) . , , 

' 1 + — ) A uu cos (^« - <>*) (2- 4 °) 



X 



V 9P J cos h kh 



Letting 




cj2 V 9P ) coshkh 



cosh k(z + h) 
u) 2 coshkh y gp J I S mhk(z + h) 




(2.41) 



(2.42) 



we get 



x' 2 z' 2 



a 2 + V =1 ^ 
The particle trajectory at any depth is an ellipse. Both horizontal (major) and vertical 

(minor) axes of the ellipse decrease monotonically in depth. The minor axis diminishes 

to zero at the seabed, hence the ellipse collapses to a horizontal line segment. In deep 

water, the major and minor axes are equal 

, gkA / Tk 2 \ kz 
a = 6 = ^-— 1 + )e kz % 2.44 

u 2 V gp J 

therefore the orbits are circles with the radius diminishing exponentially with depth. 
Also we can rewrite the trajectory as 

Tk 2 \ cosh k(z + h) . , . 

x' = 1 + V4l — " sin ( ut ~ kx o) 2 - 45 

cosh kh 





Tk 2 


U) 1 \ 


9P 


ak t (l + 


Tk 2 




gp 


' = and 


z' = 



Tk 1 \ sinh k(z + h) . . , 7i\ ,„ 

' sm{ut - kx - -) (2.46) 



cosh A; /i 2 
When cot — kx = 0, x' = and z' = b. A quarter period later, cot — k = it /2, x' = a 
and 2' = 0. Hence as time passes, the particle traces the elliptical orbit in the clockwise 
direction. 



2.5 Energy and Energy transport 



K 9 



J dz(u 2 + w 2 ^) (2.47) 



Beneath a unit length of the free surface, the time-averaged kinetic energy density is 

»o 

2 

whereas the instantaneous potential energy density is 

E P = \pgC 2 + = \ P gC 2 + T (VTTC 2 - l) = \ P gC + T£ (2.48) 

Hence the time-average is 

E P = \pge + \~Q (2-49) 
Let us rewrite (2.25) and (2.26) in (2.49): 

u = R [ 9kA (l I Tk2 ) COshk ( Z + h ) c ikA c -iut (2 50) 

\ u \ gp J cosh kh J 

w = ^ JzlgM ( x + T#) sinh *(* + h ) e *A e -™ (2 . 51) 

[ a; \ r/p / coshkh J 
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Then 



p fgkAV ( Tk 2 \ 2 smh2kh p ( gkA\ 2 f Tk 2 \ 2 sinh kh 
4 \ u) ) \ gp ) 2k cosh 2 kh 4 \ u J \ gp J k cosh k h 
pgA 2 f l + Tk 2 V gkt<mhkh _ pgA 2 / | Tlf\ ^ ^ 



4 V gp J 4 V gp 

after using the dispersion relation. On the other hand 



pgA 2 ( Tk 2 \ 
E v = ^7- 1 + (2-53) 



4 V pg 

Hence the total energy density is 

Note that the total energy is equally divided between kinetic and potential energies; this 
is called the equipartition of energy. 

We leave it as an exercise to show that the power flux (rate of energy flux) across a 
station x is 



dE_ 



[ pudz- T( X Q = -p [ $ t $ a dz - TCCt = Ec 9 (2.55) 

J-h J-h 

where c g is the speed of energy transport , or the group velocity 



_du) _c / ffi + 3 2kh 1 _ c / ^ + 3 2kh 
Cg ~dk~2\ %7I + S inh2kh j 2 | + sinh 2kh ' 

For pure gravity waves, k/k m <C 1 so that 

c s = 7T 1 + —TT^TT ( 2 - 57 ) 



2 V sinh2A;/i 

where the phase velocity is 



c= A /|tanh feft (2.58) 



In very deep water &/i ^> 1, we have 



C » = 2 = 2 Vf (2 ' 59) 
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The shorter the waves the smaller the phase and group velocities. In shallow water 
kh < 1, 

c g = c= ^/gh (2.60) 

Long waves are the fastest and no longer dispersive. 
For capillary-gravity waves with kh 1, we have 

3 



? + 3 



2 a.4.1 



C 

2 



where 



Note that e„ 



^ + 1 



<g | Tfc 3 



P 



c when k = k m , and 



_ 2tt_ /p^ 

m \ V T 1 
V J 



(2.61) 



(2.62) 



c 9 < c, 



if k<k r 



In the limit of pure capillary waves of A: ^> k m , c g 



c g = c/2 as in (2.59). 



(2.63) 

3c/2. For pure gravity waves 



3 Wave resistance of a two-dimensional obstacle 

Ref: Lecture notes on Surface Wave Hydrodynamics Theodore T.Y. WU, Calif. Inst. Tech. 

As an application of the information gathered so far, let us examine the wave resis- 
tance on a two dimensional body steadily advancing parallel to the free surface. Let the 
body speed be U from right to left and the sea depth be constant. 

Due to two-dimensionality, waves generated must have crests parallel to the axis of 
the body (y axis). After the steady state is reached, waves that keep up with the ship 
must have the phase velocity equal to the body speed. In the coordinate system fixed 
on the body, the waves are stationary. Consider first capillary -gravity waves in deep 
water A* = A/A TO = 0(1) and kh S> 1. Equating U = c we get from the normalized 
dispersion relation 

where U* = U/c m . Hence 

A* - 2c^A* + 1 = = (A* - A* a )(A* - A* 2 ) 
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which can be solved to give 



A* 2 



cl ± (ci - 1) 



1/2 



and 



A*i 



1 

A* 2 



(3.2) 



(3.3) 



Thus, as long as c* = £/* > 1 two wave trains are present: the longer gravity wave 
with length A* 1? and the shorter capillary wave with length A* 2 . Since c gi < c = U and 
c g2 > c = U, and energy must be sent from the body, the longer gravity waves must 
follow, while the shorter capillary waves stay ahead of, the body. 

Balancing the power supply by the body and the power flux in both wave trains, we 

get 

Rc= (c- c gi )E 1 + (c 92 - c)E 2 

Recalling that 



(3.4) 



we find, 



1-^ = 1-1 I 

c 2 



Cg = l Xj + 3 
c 2 A 2 + 1 

\ 1 



1/A* 1 _ 1/A* 
Xl + lJ 2 A, + 1/A* ~ 2 2c 2 



For the longer wave we replace c g /c by c 9atl /c* and A* by A*i in the preceding 
equation, and use (3.2), yielding 



"9*1 



4x1/2 



(3.5) 



Similarly we can show that 



"ff*2 



1 = (1 - c;') 



Since 



— ( 1 + X! 



we get finally 



4x1/2 



"9*1 



2 A*! y A*i 



p#A 2 A* 2 c 2 , 



= \p9 {K 2 A{ + A*^ 2 ) (4 - l) 1 / 2 = ^ (A* 2 A 2 + A,-^ 2 ) {Ut - If 



1-2 



(3.6) 



(3.7) 



(3.? 
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Note that when U* = 1, the two waves become the same; no power input from the body 
is needed to maintain the single infinite train of waves; the wave resistance vanishes. 
When U* < 1, no waves are generated; the disturbance is purely local and there is 
also no wave resistance. To get the magnitude of R one must solve the boundary value 
problem for the wave amplitudes Ai,A 2 which are affected by the size (relative to the 
wavelengths), shape and depth of submergence. 

When the speed is sufficiently high, pure gravity waves are generated behind the 
body. Power balance then requires that 

R=(l-*L)E=W*-( 1 -- kk ) (3.9) 
V UJ 2 \2 smh2kh) V ; 

The wavelength generated by the moving body is given implicitly by 



U ( tanh kh 



1/2 



(3.10) 



\fgh~ \ kh 

When U ~ y/gh the waves generated are very long, kh <C 1, c g — > c = y^gh, and the 
wave resistance drops to zero. When U \/g~h, the waves are very short, fc/i> 1, 



For intermediate speeds the dependence of wave resistance on speed is plotted in figure 
(3). 



4 Narrow- banded dispersive waves in general 

In this section let us discuss the superposition of progressive sinusoidal waves with the 
amplitudes spread over a narrow spectrum of wave numbers 

pec poo 

C{x,t)= \A(k)\cos(kx-Lut-9 A )dk = $t A(k)e ikx - iuJt dk (4.1) 
Jo Jo 

where^4(A;) is complex denotes the dimensionless amplitude spectrum of dimension 

(length) 2 . The component waves are dispersive with a general nonlinear relation uj{k). 

Let A(k) be different from zero only within a narrow band of wave numbers centered at 

k . Thus the integrand is of significance only in a small neighborhood of k . We then 
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Figure 3: Dependence of wave resistance on speed for pure gravity waves 



approximate the integral by expanding for small AA; = k — k and denote u c 
co' = co'{k ), and to" = u"(k ), 



co(k ), 



C 



ik x—iuj, 



r-oo 

ot / A(k) e ^*-^-«o)t dk 
Jo 



ik x—iu) t 



dk A{k) exp 



iAkx - i [ co' Ak + -J^Akf ) t + 



{A(x } t)e ikoX ~ iuJot } 



where 



A(x,t) 



dk A(k) exp 



iAkx - i ( u)' Ak + -u"(Ak) 2 ) t + 



(4.2) 



(4.3) 



Although the integration is formally extends from to oo, the effective range is only 
from k — (Ak) m to k + (Ak) m , i,.e., the total range is 0((Ak) m ), where (Ak) m is the 
bandwidth. Thus the total wave is almost a sinusoidal wavetrain with frequency co and 
wave number k , and amplitude A(x,t) whose local value is slowly varying in space and 
time. A(x,t) is also called the envelope. How slow is its variation? 
If we ignore terms of (AA;) 2 in the integrand, (4.3) reduces to 



A(x,t) 



dk A(k) exp [iAk(x - uj' t)\ 



(4.4) 
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Clearly A = A(x — to' t). Thus the envelope itself is a wave traveling at the speed to' 
This speed is called the group velocity, 

du) 



(4.5) 

ho 



Note that the characteristic length and time scales are (A& m ) _1 and (uj' /S.k m )~ 1 respec- 
tively, therefore much longer than those of the component waves : k~ Y and a;" 1 . In other 
words, (4.3) is adequate for the slow variation of A e in the spatial range of Ak m x = 0(1) 
and the time range of u)' Ak m t = 0(1). 

As a specific example we let the amplitude spectrum be a real constant within the 
narrow band of k — K,k + k, 

( = A I ° + e lkx ~ iu){k)t dk, k <C k (4.6) 

J k —n 

then 

/K 
fitgikoZix-Cgt) _|_ . . . 
■K 

_ 2 A sin k(x — Cgt) ^ koX _ iuJo t _ ^ e ik x-iuj t u y\ 

where £ = k — k /k and 



A = 2 A sin k(x - c g t) 

(X - Cgt) [ ' ] 



as plotted in figure (4). 

By differentiation, it can be verified that 

dA dA 



Multiplying (4.9) by A*, 

BA BA 

and adding the result to its complex conjugate, 

. dA* . dA* n 
^ + C ^ = °' 

we get 
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Figure 4: Envelope of waves with a rectangular band of wavenumbers 



We have seen that for a monochromatic wave train the energy density is proportional 
to | ^4 1 2 . Thus the time rate of change of the local energy density is balanced by the net 
flux of energy by the group velocity. 

Now let us examine the more accurate approximation (4.3). By straightforward 
differentiation, we find 



OA 

~dt 

OA 
dx 
d 2 A 
dx 2 



where 



- lL o\k )Ak- % -^^(Ak) 2 



(iAk)A(k)e ib dk 
(-(Ak) 2 ) A(k)e iS dk 



S = Akx- u' Ak t - -u"(Ak) 2 1 



A(k)e iS dk 



is the phase function. It can be easily verified that 



dA ,dA 

-m +UJ °dx~ 



iu" d 2 A 



{All) 



(4.12) 



2 dx 2 

By keeping the quadratic term in the expansion, (4.12) is now valid for a larger spatial 
range of {Ak) 2 x = 0(1). In the coordinate system moving at the group velocity, £ = 
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x — Cgt, t = t, we easily find 



dA{£,r) OA OA 8A{£,t) OA 



dt dr 9 <9£ ' dx dx 

so that (4.12) simplifies to the Schrodinger equation: 



OA iu" d 2 A 



dr 2 &e 

By manipulations similar to those leading to (4.10), we get 



(4.13) 



d\A\ 2 iu" d ( BA 8A*\ 

Thus the local energy density is not conserved over a long distance of propagation. 
Higher order effects of dispersion redistribute energy to other parts of the envelope. 
For either a wave packet whose envelope has a finite length ( A(±oo) = 0), or for a 
periodically modulated envelope (A(x) = A(x + L)), we can integrate (4.14) to give 

^-f\A\ 2 dt = (4.15) 

where the integration extends over the entire wave packet or the group period. Thus 
the total energy in the entire wave packet or in a group period is conserved. 



2 Two- Dimensional Tsunami 

Tsunamis are the water waves generated by submarine earthquakes. If the seafloor 
displacement is known in the area of the earthquake, the water-wave problem is a purely 
hydrodynamic one. Unfortunately, direct measurements near the epicenter are too dif- 
ficult to make, and a good deal of effort has been centered on using water-wave records 
measured at larger distances from the epicenter to infer roughly the nature of tectonic 
movement. Hence, there has been considerable theoretical studies on water waves due 
to a variety of ground movement. 

As an introduction, let us consider an open ocean of constant depth and assume that 
there is no wind and the tectonic disturbance on the bottom is independent of y. The 
problem is two dimensional in the x, z plane. Thus, the velocity potential <&(x,z,t) 
satisfies 
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On the free surface the following conditions hold: 

£ = , = (2.1.2a) 

— +g( = 0, z = 0, (2.1.2b) 

where the atmospheric pressure is assumed to be absent, the two conditions can be 
combined to give 

1iF + 9 a7 = ' z = < 2 - L2c > 

Let the seafloor be denoted by z = —h + H(x,t). If the ground motion is known, 
continuity of normal velocity gives 

Within the framework of linearization we assume that the amplitudes of H, dH/dt and 
dH/dx are small so that the quadratic term is negligible; hence 

d$ dH 

— = — = W(x, t) on z ^ -h . (2.4) 

^ dt v ' 7 v ' 

From (3.4) two initial conditions must be further prescribed, i.e., <&{x, 0,0) and C( :E ? 0) 
on the free surface, but nowhere else, because time derivatives appear only in the free- 
surface conditions. 

What is the physical significance of <&(x, 0, 0)? Assume that, before t = 0, all is 
calm, but at t = an impulsive pressure P a (x,t) = IS(t) is applied on the free surface. 
Integrating Bernoulli's equation from t = 0— to t = 0+, we obtain 

r 0+ i r 0+ i 

$(£,0,0+) - $(z,0,0-) + / 0Cd< = — / S(t)dt = — . 

Jo- P Jo- P 

Since <&(x,0, 0— ) = and £ must be finite, we obtain $(:r,0,0+) = //p. Thus, the 
initial value of $ represents physically an impulsive pressure acting on the free surface 
at an instant slightly earlier than t = + . 

Equations (2.1), (3.4), and (2.4) now define a boundary-value problem which formally 
resembles that of a simple harmonic case. For any finite t it is expected that no motion 
is felt at a great distance from the initial disturbance so that &(x, t) — >■ as \x\ — > oo, 
which implies that <£" — ^ as — ^ oo. Since the region does not involve any finite bodies, 
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the problem can be readily solved by applying the exponential Fourier transform with 
respect to x, defined by 



OO 1 poo 



f(k) = I e~ ikx f(x) dx , f(x) = ±- I e ikx f(k) dk . (2.5) 

J-oc llT J-oo 

By taking Fourier transform of Laplace's equation with respect to x, we get the 
general solution 

(j) = Asinh k(z + h) + Bcoshk(z + h) (2.6) 
Fourier transforms of the boundary conditions give 

$ ft + 0$* = O, z = (2.7) 

on the free surface and 

4> z = W, z = -h (2.8) 
on the mean sea bed. Let the free surface be undisturbed initially, 

o, o) = <MM, 0) = (2.9) 

We get from (2.7), 

B tt + gkA = (2.10) 

and from (2.8) 

kA cosh kh - kB sinh kh = W (2.11) 

It follows that 

B tt + uj 2 B = ^- (2.12) 
cosh/c/i 

where lo = V 'gk tanh kh, subject to the initial conditions 

£(Jfc,0)=0, and B ( (/c,0)=0 (2.13) 

For the special case where W(x,t) = Ho(x)S(t), integation across the delta function 
gives the the impulsive vertical displacement 



Wdt = H (x) (2.14) 
The initial-value problem can be replaced by 

B tt + u) 2 B = t>0+ (2.15) 
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fl(fc,0+) = 0, and B t (k,0+) = (2.16) 

cosh/c/i 

which can be readily solved fo B{t). Afterwards $(&, t) follows from, from which (,{k,t) 
is found. By invese transform, we finally get £(x,t). 

We leave it as an exercise to show that the free surface is given by 

'^C^^^ (2 - i7 » 

which can be written as 

i poo fjn(k) i r°° 

( = — dk °;/, e ik ( x cos ut = — / dk( (k)e ikx cos cot. (2.18) 

27T ./.oo COSh K/J, 27T J.oo 

where 

Clearly, cosa>£ is even in k. 

In general, we can split H (x) into even and odd parts with respect to x : Hq and 
Hq . It follows from the definition of Fourier transform that 

poo poo 

H (k) = 2 dx cos kxH^x) — 2i / dx sin kxH^ (x) , 
Jo -/o 

= H e (k) + H° (k) 

where Hq is real and even in k and Hq is imaginary and odd in k. 

For simplicity, let us consider H (x) (hence (o(x)) to be even in x. The case of odd 
Hq{x) can be dealt with similarly. Eq. (2.18) may be written 

1 f 00 

((x,t) = — dk Co cos kx cos tot 
k Jo 

i roc 

= — Re / dkCo[e i{kx - ut) + e iikx+u3t) ] . (2.20) 
2tt Jo 

The first and second terms in the brackets above represent right- and left-going waves, 
respectively. 

2.1 Asympotic behavior at large time 

For a better physical understanding, approximations are necessary. At large t we can 
employ the method of stationary phase devised by Kelvin. Heuristically, the idea is as 
follows. 
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Consider the integral 

I(t) = I fe lt9 dk (2.21) 

J a 

where / and g are smooth functions of k. When t is large, the phase tg of the sinusoidal 
part oscillates rapidly as k varies. If one plots the integrand versus k, there is very little 
net area under the curve due to cancellation unless there is a point at which the phase 
is stationary, that is, 

g'(k)=0, k = k . (2.22) 

In the neighborhood of this stationary point the oscillating factor of the integrand of 
Eq. (2.23) may be written 

ex P {it[g(k) - g(k )}} . 

The real part of exp{it[g(k) — g(ko)]} varies slowly, as sketched in Fig. 5, while the 
imaginary part slowly crosses the k axis at k = k . Therefore, a significant contribution 
to th For a better physical understanding, approximations are necessary. At large t we 
can employ the method of stationary phase devised by Kelvin. Heuristically, the idea is 
as follows. 

Consider the integral 

I(t) = ! fe lt9 dk (2.23) 

J a 

where / and g are smooth functions of k. When t is large, the phase tg of the sinusoidal 
part oscillates rapidly as k varies. If one plots the integrand versus k, there is very little 
net area under the curve due to cancellation unless there is a point at which the phase 
is stationary, that is, 

g'(k)=0, k = k . (2.24) 

In the neighborhood of this stationary point the oscillating factor of the integrand of 
Eq. (2.23) may be written 

e**™ exp{it\g(k) - g{ko)]} . 

The real part of exp{it[g(k) — g(ko)]} varies slowly, as sketched in Fig. 5, while the 
imaginary part slowly touches the k axis at k = k . Therefore, a significant contribution 
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Figure by MIT OCW. 



Figure 5: The real part of exp{it[g(k) — g(ko)]}. 

to the integral can be expected from this neighborhood. If we approximate g{k) by the 
first two terms of the Taylor expansion 

g(k) = g(k ) + ^(k-k ) 2 g"(k ), 



then the integral may be written 



dk exp 



-i(k - k )Hg"(k ) 



where the limits (a, b) have been approximated by (—00,00). Using the fact that 

'7T\l/2 



„±itk 2 



dk 



we finally have 



I ^ e U9{ko) f(k ) 



(7) 



2?r 



= ±17t/4 



1/2 



3 ±i7r/4 



f2.25l 



.*lff"(*o)L 

where the + sign is to be taken if g"(ko) < 0, and minus ifif g"(ko) > 0, . It can be 
shown by a more elaborate analysis that the error is of order 0(t~ r ). Also if there is no 
stationary point in the range (a, b), the integral is at most of order 0(t~ r ). This and 
other information can be found in Stoker (1957) or Carrier, Krook, and Pearson (1966). 

Returning to Eq. (2.20), we need certain properties of the dispersion curve as sketched 
in Fig. 6. Consider x > 0. For the first integral 



g (k) = k--u, 

it may be seen from Fig. 6(6) that there is a stationary point at 



x 

7 



co'(k ) = C g (k ) iff< {ghf 2 



(2.26) 
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Figure by MIT OCW. 

Figure 6: Variations of u, d', and to" with k. 

In the same interval (0, oo) of k, there is no stationary point for the second integral. It 
follows from Eq. (2.25) that 



c 



2ty 


1/2 


[t\u"(k )\_ 


COS 





7T 



k x - uj(k )t + — 



+ 0(r 1 ), x<{gh) 1/2 t, 
where use is made of the fact that co"(k) < [Fig. 6(c)], and 

C^Or 1 ), x>(gh) l / 2 t. 



(2.27) 



(2.28) 



Now let us examine the physics represented by Eq. (2.27). An observer moving at a 
certain speed x/t lower than {gh) 1 / 2 sees a train of sinusoidal waves of wavenumber k 
[and frequency co(ko)] whose group velocity equals x/t. The amplitude of the wavetrain 
decays as 0(t~ 1//2 ). For large x/t we see from Fig. 6(a) that ko is small, hence, a faster 
moving observer sees longer waves which are also of larger amplitude since (|a/'(&o)|) 1//2 
is less. The precise shape of H (x) affects H Q {k) hence Co(&) an d the amplitude of the 
dispersed waves. 

Summing up the views of many observers for the same t, we obtain a snapshot of the 
free surface (see Fig. 7). Thus, at a constant t, long waves are found toward the front 
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and short waves toward the rear. Now consider the snapshot at a later time t 2 > t\. Both 
observers have now moved to the right. The spatial separation, however, has increased. 
In particular, let £1 ~ £2 so that between them k,u) ~ const. The total extent of a 
monochromatic wavetrain with k,u) now stretches with increasing t, implying that wave 
crests are created in the course of propagation. 

To follow a particular wave crest at its phase speed, an observer must travel at a 
varying speed since k and C(k ) do not remain constant as the crest moves into new 
territory. However, if one moves at the group velocity of the waves of length 2ir/k , 
one only sees sine waves of this length catching up from behind and then running away 
toward the front, since their phase velocity exceeds the group velocity. 

A similar picture exists for the left-going disturbance. 

2.2 Leading waves of a tsunami due to symmetric vertical dis- 
placement 

The fastest waves correspond to k ~ and move at the speed near (gh) 1 ! 2 . In the 
neighborhood of the wave front, g'(k) ~ x/t — (gh) 1 ! 2 is small, and the phase is nearly 
stationary. Furthermore, uj"(k) ~ —(gh) l l 2 h 2 k is also very small and the approximation 
by the stationary phase method is not valid. A better approximation is needed (Kajiura, 
1963). 




Figure by MIT OCW. 



Figure 7: Space-time plot of dispersive waves between two moving observers. 
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Since k ~ 0, we expand the phase function for small k as follows: 

9(k) = k- t -u{k)^k{^)-{ghfl^ W 

(gh)V* - - 



6 



k 



~ - M 1/2 



(2.29) 



Near the leading wave, x/t — {gh) 1 / 2 can be zero; we must retain the term proportional 
to A; 3 . Again, only the first integral in Eq. (2.20) matters so that 



c 



2tt 



dk Q (k) cos(kx — ut) + O 



2 <o(0) J hf'H] - 



{ghfHH 



6 



k 3 > dk 



where use is made of the fact that £o is real. With the change of variables 

. 2[x - {ghflHf 



(gh)V 2 hn 



and k[x — (gh) l ' 2 t] = Za . 



the integral above becomes 



where 



(2) 1/3 Co e (0) 
2n({gh) 1 / 2 hH) 1 / 3 



C e (0) = ^ e (0) 



da cos ( Za 



a' 



Hl(x)dx 



(2.30) 



(2.31) 



is proportional to the total area of the seafloor displacement. The integral in (2.30) 
above can be expressed in terms of Airy's function of Z: 



1 f°° ( a 

Ai{Z) = — I da cos Za + — 

7T 7n V 3 



(2.32) 



Thus, we have 



(ghy/ 2 h 2 t 



1/3 



;i?o(0)Ai 



{ghy/ 2 h 2 t 



1/3 



[x - (ghflh] 



(2.33) 



Ai(Z) is oscillatory for Z < and decays exponentially for Z > 0. Its variation is shown 
in Fig. 8. 
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Figure by MIT OCW. 

Figure 8: Leading wave due to a symmetrical surface hump or trough. The ordinate is 
C[(gh)^hH/2]^[Qm-\ see Eq. (2.33). 

The physical picture is as follows: For a fixed t, Z is proportional to x — (gh^'H 
which is the distance from the wave front x = (gh) l ' 2 t. At a fixed instant the amplitude 
is small ahead of the front, and the highest peak is at some distance behind. Toward 
the rear, the amplitude and the wavelength decrease. Since Z is proportional to t~ 1//3 , 
the snapshots at different times are of the same form except that the spatial scale is 
proportional to the factor t 1 ' 3 , meaning that the same wave form is being stretched 
out in time. During the evolution the amplitude decays as t -1 ' 3 while the rest of the 
wavetrain decays as t -1 ' 2 . Thus, the head lives longer than the rest of the body. Note 
that the amplitude of the leading wave is proportional to Hq(0) which is equal to the 
total area of the initial displacement Hq(x). 

2.3 Tsunami Due to Tilting of the Bottom 

Among the many features of tsunamis as recorded near a coast, two have been frequently 
(but not always) reported (Shepard, 1963). One feature is that the arrival of a tsunami 
is often preceded by the withdrawal of water from the beaches, and the other is that 
the first crest may not be the largest. In this section we shall show an idealized model 
which reproduces these features qualitatively. 

Again any H (x) can be thought of as the sum of Hq(x) and Hq(x) which are odd and 
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even in x, respectively. By linearity, the two parts may be treated separately first and 
their results superimposed later. It is easily shown that the even part Hq(x) has effects 
very similar to the previous example of symmetrical initial displacement on the free 
surface, the only difference being the factor (coshft/i) -1 which cuts down the influence 
of the short waves. We shall, therefore, only focus our attention to the odd part. 
Let us introduce 

H° (x) = ^ (2.34) 
so that H$(k) = ikB(k). Since Hftik) is odd, B must be real and even in k; hence, 

r,ikx -y 



c 



dk 



1 

27 

1 d 

2tt dx 



27T dx 



-ikB(k)^-(e luJt 
cosh kh v 'T 



-iu)t 



dk- 



cosh kh 
dk- 



1 

"2 { 



B(k)^-(e iuJt + e~ iuJt ) 



-B(k)(e 



iust 



-iuit 



)■ 



(2.35) 



cosh kh 

For large t and away from the leading waves, the integrals can be dealt with by the 
stationary phase method just as before, and many of the same qualitative features should 
be expected. Let us only look at the neighborhood of the leading waves propagating 
to x > 0. Again, the second integral dominates and the important contribution comes 
from the neighborhood of k ~ 0. Hence 



Re 



cosh 
Re 5(0) 



B(k) 



dk e ikx e- iut 



JO 



ReB(0) j dk exp ^ ^k[x - (gh) 1/2 t] + ^(gh) 1/2 h 2 kh 



ttB(O) 



nl/3 



_(ghy/ 2 h 2 t_ 



Ai 



nl/3 



_( g hy/ 2 hn_ 



[x - (ghflH] 
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as discussed earlier. Differentiating with respect to x, we have 



B(0) 
2 

B(0) 



{gh) l l 2 h 2 t 
2 

(gh)V 2 hH 



1/3 



dx 



Ai 



{ghy/ 2 h 2 t 



1/3 



[x - (ghflH] 



2/3 



Ai' 



(gh)V 2 hH 



1/3 



[* - (ghf'h] 



(2.36) 



where 



Ai^)^Ai(Z). 



The leading wave attenuates with time as £~ 2 / 3 which is much faster than the case 
of a pure rise or fall (where ( ~ t~ 1//3 ). This result is due to the fact that the ground 
movement is half positive and half negative, thereby reducing the net effect. The function 
Ai'(Z) behaves as shown in Fig. 9. Note that 



B(0) 



B{x) dx 



dx / H°(x')dx' 



xHq(x) dx . 



Thus, if the ground tilts down on the right and up on the left, B(0) > and the wave 
front propagating to the right is led by depression of water surface (hence withdrawal 
from a beach). The subsequent crests increase in amplitude. On the left side, x < 0, the 
wave front has the opposite phase and is led by a crest. If, however, the ground tilt is 
opposite in direction, that is, down on the left and up on the right, then the right-going 
wave front should be led by an elevation. 

Kajiura pointed out that retaining the terms gk 3 h 2 in co{k) implies keeping dispersion 
to the lowest order, and the same results, Eqs. (2.33) and (2.36), may be obtained 
alternatively by invoking the long- wave approximation at the outset, which is clearly 
appropriate far away from the source. It will be shown in Chapter Twelve that such 
an approximation is given by the linearized Boussinesq equations which are, in one 
dimension, equivalent to 

r)' 2 r f FPC h' 2 fiV\ 

(2.37) 



d 2 C , fd 2 C h 2 d\ 
gh 



dt 2 



dx 2 



3 dx 4 
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Figure by MIT OCW. 

Figure 9: Leading wave due to antisymmetric ground tilt ([B(0)]~ 1 [(gh) 1 / 2 h 2 t/2] 2 / 3 , see 
Eq. (2.36). 

3 Radiation of surface waves forced by an oscillating 
pressure 

We demonstrate the reasoning which is typical in many similar radiation problems. 
The governing equations are 



V (j> = (f) xx + 4>yy =0, — OO < Z < 0. 



(3-1] 



with the kinematic boundary condition 



Ct, z = 



(3.2) 



and the dynamic boundary condition 



P 



9 



(3.3) 



where p a is the prescribed air pressure. Eliminating the free surface displacement we 
get 



{Pa 



9 



z = 0. 



(3.4) 



Let us consider only sinusoidal time dependence: 



p a = P(x)e 



-ILOt 



(3.5) 
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and assume 

(f){x, z, t) = $(rc, z)e~ iut , ({x, t) = r ? (x)e^ i 
then the governing equations become 

V 2 $ = + $ yy = 0, - oo < z < 0. 

= —iujrj, z = 



and 



$ 2 $ = — P(z), z = 0. 

9 99 

Define the Fourier transform and its inverse by 



/(<*)=/ dxe-^f(x), f(x) = ^l dae iax f{a) 

J —oo ^ J -oo 

We then get the transforms of (3.1) and (3.4) 

$ M - a 2 4> = 0, 2 < 

subject to 

The solution finite at 2 ~ — oo for all a is 

$ = 4 e M* 

To satisfy the free surface condition 

, 4 oj 2 A iuP(a)P(a) 
\a\A A = 



hence 



.4 



P# 

i£ii-P(a) 
PS 



\a\ — co 2 /g 



or 



2tt ./.^ |a| - w 2 /# 

/ dae iax e^ z / dx' e~ iax ' Pix')—, ^— , 

P# 2tt J.^ J.^ \a\ - u) 2 /g 

— dx'P(x')— (°° dae^-^e^ , , 1 „ , 
P£ J-co 2^ ./.^ |a| - w 2 /^ 
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Let 



we can rewrite (3.13) as 



.J 2 



k = — (3.14) 
9 



$ = — / dx P(x)- I dae az — — - — — (3.15) 

99 J-oo ^ Jo oi-k 



The final formal solution is 



If we chose 



"V" f dx'Pix') 1 - rdae ^ ^-^ (3.16) 
P9 J-oo k Jo a-k 



P(x') = P 8{x') (3.17) 



then 



iuPn 1 f°° , _, r cos(m;) 

§^g(x,z ) = / dae az 3.18 

P9 7T J a - k 

is clearly the response to a concentrated surface pressure and the response to a pressure 
distribution (3.16) can be written as a superposition of concentrated loads over the free 
surface, 

/oo 
dx' P(x')g(x-x',z). (3.19) 
-oo 



iuPo^t 1 r , „az C0S ( aX ) 



where 

g(x,z,t) = ^^e~ mt - I dae at % ~"^ (3.20) 
pg 7T Jo OL-k 

In these results, e.g., (3.20), the Fourier integral is so far undefined since the integrand 
has a real pole at a = k which is on the path of integration. To make it mathematically 
defined we can chose the principal value, deform the contour from below or from above 
the pole as shown in figure (3). This indefiniteness is due to the assumption of quasi 

a=k Re a Re a a=k Re a 

I \zj ^ 



a k 



Figure 10: Possible paths of integration 

steady state where the influence of the initial condition is no longer traceable. We must 
now impose the radiation condition that waves must be outgoing as x — > oo. This 
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ilm a 




► Re a 



Figure 11: Closed contour in the upper half plane 

condition can only be satisfied if we deform the contour from below. Denoting this 
contour by T, we now manipulate the integral to exhibit the behavior at infinity, and to 
verify the choice of path. For simplicity we focus attention on Q. Due to symmetry, it 
suffices to consider x > 0. Rewriting, 

99 ^ 



pg 2tt 



r 



+ 



a — k a — k 



(3.21) 



Consider the first integral in (3.21). In order that the first integral converges for 
large \a\, we close the contour by a large circular arc in the upper half plane, as shown 
in figure (11), where Qa > along the arc. The term 



is exponentially small for positive x. Similarly, for the second integral we must chose 
the contour by a large circular arc in the lower half plane as shown in figure (12). 
Back to the first integral in (3.21) 



7 a = / da (3.22) 

a — k 



r 



The contour integral is 



da — = / da — h / da — h / da — 

a — k J T a — k J c a — k J ioo a — k 

rO ^ictx gaz 

= h + + / da — 

Jioo a — K 
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Figure 12: Closed contour in the lower half plane 

The contribution by the circular arc C vanishes by Jordan's lemma. The left hand side 
is 

LHS = 2irie ikx e kz (3.23) 

by Cauchy's residue theorem. By the change of variable a = i(3, the right hand side 
becomes 

RHS = h + i I d(3- 



e -f3x e i/3z 



i/3-k 



oo e -Px e ifiz 



Hence 

h = 2me ikx e kz + i \ d(3 ~ :0 ~— (3.24) 

Now consider J 2 



i(3 - k 



and the contour integral along the contour closed in the lower half plane, 



da — = Io + + / da 



^—tax^az 



a — k Jq a — k 

Again no contribution comes from the circular arc C. Now the pole is outside the 
contour hence LHS = 0. Let a = —i/3 in the last integral we get 



00 e ~Px e -iPy 



JO 



Adding the results (3.24) and (3.26). 

h + h = 2irie lkx e kz + / d(3 



Jo 



i{3 — k —i{3 — k 



oo e -j3x 



2me ikx e kz + 2 I df3-^ — — (B cos 0y + k sin 0y) (3.27) 



(3 2 + tf 
34 



Finally, the total potential is, on the side of x > 0, 



pg \ 2ni 

= e -tot) e ikx e kz + _ d a (/Jcoa/fc + jfesin/fc) \ (3.28) 

PQ I n J f3 2 + k 2 J 

The first term gives an outgoing waves. For a concentrated load with amplitude 
P , the displacement amplitude is P /pg- The integral above represent local effects 
important only near the applied pressure. If the concentrated load is at x = x', one 
simply replaces x by x — x' everywhere. 

Note on Rayleigh's fictitious damping: If frictional dissipation is accounted 
for, there should be no wave at infinity. A way to derive the solution for dissipationless 
theory is to allow a small damping at first, obtain the solution subject to the condition 
that there is no radiation at infinity, and then take the limit of zero damping. In some 
problem it is easy to include the real damping mechancism. in some others a fictitious 
damping can be used. 

For linearized water waves we can introduce a fictitious damping term in the mo- 
mentum equation: 

-^ = —Vp-ge z -eu (3.29) 

where e is the Rayleigh damping coefficient. Introducing the velocity potential u = V$ 
we get after integration with respect to space the fictitious Bernoulli equation 

— + £+gz-e$ = (3.30) 
at p 

On the free surface it gives 

^ + ^+K-e$ = (3.31) 
ot p 

Combined with the kinematic boundary coondition we get on the free surface 

a 2 $ 9 $ a $ i dPa 

w +9 di- e ^ = -p^r (3 - 32) 

With the new term the Fourier integral has two simple poles slighlty off the real axis. 
One is in the first quadrant and the other in the third. When the limit of zero damping 
is taken a the end. they give the same residues as before, by Cauchy' theorem. 

This devise is of course artificial since in real fluid dissipation by viscosity is of a 
different form. 

We illustrate the use of Rayleigh's damping for two ship-wave problems. 
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4 Linearized equations in the moving coordinate sys- 
tem 



For simplicity let us assume that the disturbance travels at constant speed U from right 
(x' rsj oo) to left (x 1 ~ — oo) in the stationary coordinate system (x',y,z,t) where water 
at infinity is at rest. Let (x,y,z,t) be the coordinate system fixed on the disturbance, 
then the two systems are related by 



x = x ' + Ut' } y = y' } z = z', t = t'. (4.1) 



Using the chain rule 



dF(x(x', t), y, z, t) _ dF dt dF dx _ dF dF 

dt' ~~dtdt' + ~dxdt'~~dt + ~dx' 

dF(x{x',t) i y,z,t) _ dF dx _ dF 

dx' dx dx' dx 

dF _dF dy _ dF dF _ dF dz _ dF 

dy' dy dy' dy ' dz' dy dz' dz 

The dynamic and kinematic boundary conditions on the free surface z = become, 

respectively: 

Qt + UQ. + gC + & = (4.2) 

P 

where e is the Rayleigh damping factor, and 

Q + UQ = (4.3) 

Inside the fluid we still have 

®xx + ®yy + <&zz = 0, Z < 0. (4.4) 

Assuming steady state, (4.2) and (4.3) reduce to 

U* a +gC + e4> = 2 = 0, (4.5) 

P 

and 

UQ = z = 0, (4.6) 

which can be combined to 

+ = * = (4.7) 

9 9 P9 dx 

In the following examples we restrict to deep water so that f -)■ as z ^ -oo. 
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5 Two-dimensional waves due to a traveling surface 
presssure 



5.1 Solution by Fourier transform 

The governing equations for <&(x,z) are 

$** + ®zz = 0, - h< z < 

3 ^ e * 1 <9Pa 



pU dx 
$ -> 0, 2 -)■ -oo, 



^ = 0; 



(5.1) 

(5.2) 
(5.3) 



Applying the exponential Fourier transform with respect to x, defined with its inversion 

by 



$(M) 



2n 



e lkx $(x,z)dx, <&(x,z) = -^- I e lkx $(k,z)dk 



we get a two-point boundary value problem for $ 



r$ = 0, 



dz 2 

As an example we take 



/i < z < 



p a (k), z = 0; 



p a {x) = P e x2 / 4L \ so that p a (k) 



P L k 2 L 2 



The solution to (5.5) is 



$ = Ae \k\z 



(5.4) 

(5.5) 
(5.6) 

(5.7) 
(5.8) 



which vanishes at great depth. To satisfy the surface boundary condition we must have, 

-ikp a 



3 |7i , 2 , * eA: 



T 



pf/ 



hence 



and 



A 



j-\h\ — hi 4- isk 



Hk,y,z) 



JL\k\ — k 2 -A- M 



(5.9) 



(5.10) 



(5.11) 



u 2 1 
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By inverse transform the potential is 



i r 00 — — p (k)e^ z 
*(x, y,z) = - dk e ik * ■ * a \ ' (5.12) 



9tt JL I h I — £-2 i iek 

J — 00 JJ2 1^1 "> ^ (J 



The Fourier integral can be split into two: 



1 f°° 

®{x,y,z) = —J dk 




't\z 


9 \h\ Jc 2 4- 


iek 
U 


1 f° 

+ ^L dk 




k\z 


3_\h\ Z~2 _|_ 


iek 
U 


\k = k. In the second \k\ - 


= —k, hence 





1 r ite ~$Pa(k)e kz 



1 f°° 

®(x,y,z) = — dk e lkx - 

Z7r Jo 77: 



jj 2 A, ^ [/ 



1 r „ -§Pa(k)e- kz 
g_ 



Let us replace A; by -k in the second integral so that 

1 /"°° 

®{x,y,z) = — dk e 
2tt Jo 



1 f 00 „ -%Pa{k)e kz 
thk-k* + f 
_ lkx #Pa(k)e kz 
&k-k 2 -f 



+ — / dk e 



2tt P U J k-£-% 



(5.13) 



dk e lkx — pu ^ ' -r (5.14) 



27rp(7 Jo - £ + f 

Use has been made of the fact p a (k) given by (5.7) is even in k. 

5.2 Asymptotic solution in the far wake 

For brevity we shall write (5.15) as 

®(x,y,z) = $i + $ 2 (5.16) 

where $i and $ 2 represent the first and second integrals respectively. Let us examine 
each integral in turn in the complex k plane. The first integrand has a pole in the 
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Figure by MIT OCW. 

Figure 13: Contours in the complex k plane: (a) for $! and x < 0.,(b) for $ 2 and x < 0, 
(c) for $! and x > 0, (d) for $ 2 and x > 0. 



first quadrant; the second has one in the fourth quadrant. If no Rayleigh damping 
were assumed, both poles would be on the path of integration, rendering the integrals 
undefined, unless the radiation condition is added, see §3. 

Consider the first x < i.e., ahead of the ship. We replace the integral path for $ x 
by the closed clockwise contour in the fourth quadrant, as shown in Figure 5.2-a. The 
contour consists of the orginal path, hence $i, an integral along a circular arc C_ with 
very large radius, and the positive imaginary axis, all with the same integrand, 

Pa(k)e kz A , % f f , f° \ Jl Jkx p a (k)e kz 



dk e ikx — 
2npU J k 



9 



Mr 

u 



$1 



2?rpU 



dke 



9 

u 2 



it 
u 



(5.17) 



Note that the contour is closed in the fourth quadrant where ki < to ensure convergence 
for x < 0, since 

(5 18) 



^ik r x g k{X 



Let us show first that the integral along the imaginary axis has only local effects, 



dke 



ikx 



Pa(k)e 



k~ 



h 9 «e 

^ U 2 U 



(5.19) 
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Let k = —iK where Kb is real and positive, we get 

v^C'ftH*' (5 20) 

This integral clearly diminishes to zero as x — >■ oo. Next, the integral along the arc C_ 
vanishes by Jordan's lemma. Since no pole is inside the contour, the contour integral 
is zero by Cauchy's residue theorem. Hence $1 has only local effect and diminishes to 
zero as far ahead of the moving pressure. 

Similarly for $ 2 we choose a closed contour in the first quadrant where ki > as in 
Figure 5.2-b, since 

e -ifcr _ e -i* r * e *<* (5.21) 

Again there is no pole inside the contour, hence the contour integral vanishes. Since 
nothing comes from the circular arc, and the integral along the imaginary axis has 
only local effect, $ 2 has only local effects also, In summary, far ahead of the moving 
disturance, wave motion ($ = $1 + $ 2 ) is negligible. 

Let us now consider x > 0, i.e., the wake. For $i the closed contour must now be in 
the first quadrant where there is the pole at 

+ £ (5-22) 

see Figure 5.2-c. By Cauchy's residue theorem the (counter-clockwise) contour integral 
is 

jdk e ik \ M T^ k = -^ l9X ' U ' 2 P« (f) e9Z/U ' 2 ( 5 - 23 ) 
By repeating earlier reasoning, integrals along the arc and along the positive imarginary 
axis do not matter in the far field, hence 



2irpU 



pU Vtt 

On the other hand, for <&i the closed contour must now be in the fourth quadrant where 
there is the pole at 

see Figure 5.2-d. By Cauchy's residue theorem the contour integral in the clockwise 
direction is 

j dk e-** ] ^ 1£ = ~e--/"-p a (H) (5.26) 



2ixpU 



u 2 ' u 
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Again, integrals along the arc and along the positive imarginary axis do not matter in 
the far field, hence 

$ 2 « _ _Lf^ e -i9x/ui e -( 9 L/u>r e9Z /ui (527) 

pU y/TT 

Finally in the far wake, 

$ = $1 + $ 2 « -\- T ^= cos(gx/U 2 )e- i9L/u2)2 e 9z/u2 (5.28) 

pU V7T 

Thus one sees a train of sinusoidal surface waves of length U 2 jg following the distur- 
bance. The wave amplitude is 

P e- (9L/c/2)2 (5.29) 

which is large if the pressure variation is sharp, and small if it is flat. 
We leave it as an exercise to derive the wave drag. 

6 The physics of Kelvin's ship- wave pattern 

Refs: Explansions are due to LIGHTHILL, First ONR SYMPOSIUM on Naval Hydro- 
dynamics (1959?). 

Material here is borrowed from the lecture notes by T. Y. Wu, Caltech. 

Mathematical details can be found in Stoker: Water Waves 1957. 



The action of the ship's propeller 
Has a thrust pattern 

To which the ship reacts by moving forward, 

Which also results secondarily, 

In the ship's bow elevated waves, 

And its depressed transverse stern wave, 

Which wave disturbances of the water 

Are separate from the propeller's thrust waves. 
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R.Buckminster Fuller, Intuition- Metaphysical Mosaic. 1972. 

Anyone flying over a moving ship must be intrigued by the beautiful pattern in the 
ship's wake. The theory behind it was first completed by Lord Kelvin, who invented 
the method of stationary phase for the task. Here we shall give a physical/geometrical 
derivation of the key results 

Consider first two coordinate systems. The first r = (x, y, z) moves with ship at the 
uniform horizontal velocity U. The second r' = (x 1 , y', z) is fixed on earth so that water 
is stationary while the ship passes by at the velocity U. The two systems are related by 
the Galilean transformation, 

r' = r + Ut (6.1) 
A train of simple harmonic progressive wave 

C = 5R{Aexp[z(k-r'-cjt)]} (6.2) 

in the moving coordinates should be expressed as 

C = M{Aexp[ik-(r-Vt) -iujt]} = M{Aexp[ik-r-i(oj-k.V)t]} 

= M{Aexp\ik-r-iat}} (6.3) 

in the stationary coordinates. Therefore the apparent frequency in the moving coordi- 
nates is 

cr = u) - k • U (6.4) 

The last result is essentially the famous Doppler's effect. To a stationary observer, the 
whistle from an approaching train has an increasingly high pitch, while that from a 
leaving train has a decreasing pitch. 

If a ship moves in very deep water at the constant speed — U in stationary water, 
then relative to the ship, water appears to be washed downstream at the velocity U. 
A stationary wave pattern is formed in the wake. Once disturbed by the passing ship, 
a fluid parcel on the ship's path radiates waves in all directions and at all frequencies. 
Wave of frequency uj spreads out radially at the phase speed of c = g/u) according to 
the dispersion relation. Only those parts of the waves that are stationary relative to the 
ship will form the ship wake, and they must satisfy the condition 

a = 0, (6.5) 
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Figure 14: Waves radiated from disturbed fluid parcel 

i.e., 

to k 

a> = k • U, or c = — = — • U (6.6) 
k k 

Referring to figure 14, let O, (x = 0) represents the point ship in the ship-bound 
coordinates. The current is in the positive x direction. Any point X\ is occupied by 
a fluid parcel Q\ which was disturbed directly by the passing ship at time t\ = X\/U 
earlier. This disurbed parcel radiates waves of all frequencies radially. The phase of 
wave at the frequency to reaches the circle of radius cti where c=g/u by the deep water 
dispersion relation. Along the entire circle however only the point that satisfies (6.6) 
can contribute to the stationary wave pattern, as marked by P. Since OQi = X\ = Utx, 
Q\P = ct\ and OP = Uti • k/k, where k is in the direction of Q\P. It follows that 
AOPQi is a right triangle, and P lies on a semi circle with diameter OQ\. Accounting 
for the radiated waves of all frequencies, hence all c, every point on the semi circle 
can be a part of the stationary wave phase formed by signals emitted from Q\. Now 
this argument must be rectified because wave energy only travels at the group velocity 
which is just half of the phase velocity in deep water. Therefore stationary crests due 
to signals from Qi can only lie on the semi-circle with the diameter 0\Q\ = OQi/2. 
Thus Pi instead of P is one of the points forming a stationary crest in the ship's wake, 
as shown in figure 14. 

Any other fluid parcel Q 2 at x 2 must have been disturbed by the passing ship at time 
h = X2/U earlier. Its radiated signals contribute to the stationary wave pattern only 
along the semi circle with diameter O2Q2 = OQ2/2. Combining the effects of all fluid 
parcels along the +x axis, stationary wave pattern must be confined inside the wedge 
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O Oi Oi Qi 



Figure 15: Wedge angle of the ship wake 




o o' Q 



Figure 16: Geometrical relation to find Points of dependence 

which envelopes all these semi circles. The half apex angle (3 of the wedge, which defines 
the wake, is given by 

si ^° = mji = m (6 ' 7) 

hence (3 = sin -1 1/3 = 19.5°, see figure 15. 

Now any point P inside the wedge is on two semicircles tangent to the boundary 
of the wedge, i.e., there are two segments of the wave crests intersecting at P: one 
perpendicular to PQi and one to PQ2, as shown in figure 15. 

Another way of picturing this is to examine an interior ray from the ship. In figure 
(16), draw a semi circle with the diameter O'Q = OQ/2, then at the two intersections Pi 
and P 2 with the ray are the two segments of the stationary wave crests. In other words, 
signals originated from Q contribute to the stationary wave pattern only at the two 
points Pi and P2, as shown in figure 16. Point Q can be called the point of dependence 
for points Pi and P 2 on the crests. 

For any interior point P there is a graphical way of finding the two points of depen- 
dence Qi and Q?. Referring to figure 16, AO'QPi and AO'QP2 are both right triangles. 
Draw 0\M\ || QP\ and 2 M 2 \\ QP2 where M\ and M 2 lie on the ray inclined at the an- 
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Figure 17: Points of dependence 

gle p. it is clear that OM x = OP L /2 and OM 2 = OP 2 /2, and AM L 0'I\ and AM 2 0'P 2 
are both right triangles. Hence O' lies on two semi circles with diameters M X P\ and 
M 2 P 2 . 

We now reverse the process, as shown in figure 17. For any point P on an interior 
ray, let us mark the mid point M of OP and draw a semi circle with diameter MP. 
The semi circle intersects the x axis at two points Si and S 2 . We then draw from P two 
lines parallel to MS\ and MS 2 , the two points of intersection Qi an Q 2 on the x axis 
are just the two points of dependence. 

Let LPQ x O = IMSxO = 9 X and LPQ 2 = LMS 2 = 6 2 . then 

PS PS- 

hence 

tan 6i + tan j3 



2 tan 9j 



1 — tan 9i tan /5 
which is a quadratic equation for 9i, with two solutions: 



tan 9i I 1 ± ^/i - 8 tan 2 /? 



(6.8) 



tan 9 2 I 4 tan/? 

They are real and distinct if 

l-8tan 2 /3>0 (6.9) 

These two angles define the local stationary wave crests crossing P, and they must 
be perpendicular to PQ\ and PQ 2 . There are no solutions if 1 — 8tan 2 /3 < 0, which 
corresponds to sin (3 > 1/3 or j3 > 19.5°, i.e., outside the wake. At the boundary of the 



wake, (3 = 19, 5° and tan f3 = \/l/8, the two angles are equal 



/2 

9 X = 9 2 = tan" 1 -y- = 35°16'. (6.10) 
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Figure by MIT OCW. 



Figure 18: Diverging and transverse waves in a ship wake 

By connecting these segments at all points in the wedge, one finds two systems of wave 
crests, the diverging waves and the transverse waves, as shown in figure div-trans. 
A beautiful photograph is shown in Figure 19 

Knowing that waves are confined in a wedge, we can estimate the behavior of the 
wave amplitude by balancing in order of magnitude work done by the wave drag R and 
the steady rate of energy flux 

RU = (Ec g )r ~ {\A\ 2 c g )r (6.11) 

hence 

A ~ r 1 ' 2 (6.12) 
This estimate is valid throughout the wedge except near the outer boundaries, where 

A ~ r" 1/3 (6.13) 

by a more refined analysis (Stoker, 1957, or Wehausen & Laitone, 1960). 
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"Photograph removed due to copyright restrictions." 



Figure 19: Ships in a straight course. From Stoker, 1957.p. 280. 
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7 Three-dimensional analysis of Kelvin's ship wave 

References: J. N. Newman: Marine Hydrodynamics 
J. J. Stoker: Water Waves, 

J. V. Wehausen and E. V. Laitone Surface Waves in Handbuch der Physik: Band IX. 
Springer. 

T Y. Wu Lecture Notes on Water Waves, Calif. Inst. Tech. 

7.1 Solution by Fourier transform 

Define the double Fourier transform: 

p poo 

$(a,(3,z)= e- iax - i/3y $(x,y,z)dxdy (7.14) 

J J — 00 

and its inverse: 

®(x,y,z) = ff e iax+lPy $(a,/3,z)dad(3 (7.15) 

^ " J J — oo 

Applying the Fourier Transform to Laplace's equation and the boundary conditions we 

get 

$ 22 - = 0, where k = y V + /? 2 , z < (7.16) 

®z - (-— it®) $ = -ict—p a (7.17) 

V 9 J P9 



The solution is easily found to be 



U iap a e kz 

> — 

Formally the inverse transform is 



pg - i ea - k 



®(x, y,z) = -^ fT dad? i^- J^^ , (7-19) 



^JJ-oo pg^f-iea-k 
The Fourier transform of the kinematic free surface condition is 

iaU((a, (3) = ® z (a, (3, 0) (7.20) 

Hence the displacement is 

«*• v) = ^ /£ ^ e "' + * c^ (7 - 2i) 
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We shall take a Gaussian pressure distribution for illustration, 

Pa = Poe^ 2+ ^ 4L2 (7.22) 

The Fourier transform is 

Pa = ^V- 2 (7.23) 

7T 

From here on we shall discuss the free surface displacment only, 



C(x, y) = -2- — / / dad(3 e^y^—, (7.24) 

7.2 Asymptotic analysis of the wake 

If we let 

x = rcos6, y = rsin#; a = k cos ip, f3 = k sin ip (7.25) 

then 

ax + j3y = kr cos(# — ip) (7.26) 

(7.24) becomes 

p /" 2 1 r 27T r°° kp- k2L2 



P L 2 5 1 /** <ty r^ e ^co S (^) fce "" 2L2 (7.27) 



pgrr U 2 4tt 2 J q cos 2 ^ 7 _ _ ieg 

We now evaluate the integral in the complex plane of k = \k\e ia . The integrand has 
a simple pole at 

k = 9 + %tg (7.28) 
t/ 2 cos 2 *0 [7 2 cos ^ 

which is slightly above the real axis. Let us use Cauchy residue theorem by choosing a 
closed contour in the complex k plane. Since 

ikr cos(-0 — 9) = (i\k\r cos a — \ k\r sin a) cos(ip — 6) (7.29) 

For < [?/> — 0[<7r/2we must choose < a < tt/2 for the Fourier integral to converge. 
Thus we close the counter-wise contour in the first quadrant by adding a large circular 
arc connecting the ends of the real and imaginary axes k ~ oo and k ~ «oo, and then 
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the positive imaginary k axis, as shown in Figure 5.2. The contour integral in k is equal 
to the residue at the simple pole by Cauchy's theorem, 



+ / + 

c 



dk 



2nig 



U 2 cos 4 ip 



exp 



( sL Y 




\ U 2 cos 2 ip J 


exp 



. gr cos(ip — 9) 
U 2 cos 2 ip 



(7.30) 



The limit of e = has been taken after calculating the residue. It can be shown that 
the line integral along the imaginary axis has only local effects and dies out quickly for 
large r. The line integral along the circular arc also vanishes in the limit of large \k\ 
by Jordan's lemma. Hence the residue is the only important term representing the line 
integral along the real k axis. Subsitituting this result in (7.27), we get 

\ 2" 



((x,z) 



2tt 2 pg U 2 



dip 
cos 4 ip 



exp 



gL 



U 2 cos 2 ip 



exp 



. gr cos(V> — 9) 
U 2 cos 2 ip 



(7.31) 



For 7r/2 < \ip — 9\ < ir we must choose — n/2 < a < for the Fourier integral to 
converge. In other words, we take a clock-wise countour closed by a circular arc in the 
fourth quadrant. The pole is now outside the contour and the residue is zero. Hence 
the k integral is negligible. 

In the far wake fj^ ^> 1, we apply the method of stationary phase to (7.31). Let the 
phase function be 



cos('0 — 9) 



Its derivative is 



F^,9) 



dF_ 

dip 



cos 2 ip 



sin(-0 — 9) sin ip cos(-0 — 9) 
i 2 



(7.32) 



cos 2 ip 
cos{ip — 9) 



cos 3 ip 
tan(-0 - 9) + 2 tan ip) 



cos 2 ip 

= F [- tan(ip - 9) + 2 tan ip) 
The points of stationary phase are the roots of 

-tan(^ - 9) +2 tan?/; = 



(7.33) 



(7.34) 



or 



2 tan 9 tan 2 ip + tan ip + tan 9 = 



(7.35) 
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Note that tan# = y/x represents the ray where an observer sits. Because of symmetry 
with respect to the x axis we need to consider x > 0, y > only. There are two negative 
solutions 



tan ^ \ = -l T yi-8(^) 2 
tan ^ 2 J 

if (x,y) is inside the half-wedge tan# = y/x < l/v8, i.e., < 19.3°. Outside the wedge 
no stationary phase point exists. 
By definition (7.26), 

tanip = - (7.37) 
a 

gives the direction of the local wavenumber vector. For a point along the edge of the 
wake, y/x = l/v8, 

tan^i = tan-02 = y=, or, tpx = ip 2 = — 35°16'. (7.38) 

For a point along the ray y/x ~ 0, 

tan^ fa — ~ — oo, (7.39) 



hence i/^ ~ — |, and 



1 - a/1 - 8(y/xy 4(y/xY y 

tan^ 2 = A i ~ = 7 -40 

Ay fx 4y/x x 

hence -02 ~ —0. Wave crests corresponding to tpi are locally parallel to the the ship's 
path (the x axis). These are the diverging waves. On the other hand, wave crests 
corresponding to ip 2 are locally perpendicular to the ship's path; they are the transverse 
waves. The ranges of tan tp are 

Diverging wave : -= > tan -01 > — oo, 

v2 

Transverse wave : > tan ip 2 < 7= (7-41) 

v2 

To apply the formula of the stationary phase method we calculate F^(tp) from (7.33) 
at the stationary points: 

F^{ip m ) = F^[- ta,n(ip m - 9) + 2tan^ m ] + F[- sec 2 (ip m - 6) + 2 sec 2 -0m], m = 1,2. 
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At the stationary points the first term vanishes. We rewrite the second term as 



*W(^m) = F t 1 - tan 2 (^ m - 9) + 2tan 2 m ] = F [l — 2 tan 2 m ] 



(7.42) 



Use has been made of (7.34). It is clear from (7.41) that F^(ipi) < for the diverging 
wave and F^{ip2) > for the transverse wave. It follows that the free surface elevation 
in the wake is dominated by 



2irU 2 



i PqL 2 9 
2tt 2 pg U 2 | \gr\F^(>4>2 



1/2 exp 



( 9L V 

V U 2 cos 2 ipi J 



cos 4 ipi 



exp 



. gr cos('0i — 9) 
U 2 cos 2 ipi 



1.71 

T 



+ 



2nU 2 



gr\F ff {ip 2 )\ 



1/2 exp 



( 9L V 

y £/ 2 cos 2 ^i / 



COS 4 '02 



exp 



. gr cos(-02 — 9) m 
U 2 cos 2 "02 4 



(7.43) 



Thus the two wave systems are out of phase bu tt/2. Each wave decays with distance 
as 1/ y/r. The wave length of each wave system is 

2t\\J 2 cos ip m 



Xr 



9 



(7.44) 



For the diverging waves, cos ipi « along the centerline, where the wave length is the 
shortest and increases outwards. 

Note that near the edge of the wake, F^ = 0. The present approximation breaks 
down. A better approximation can be made (see Lamb, Stoker). 
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8 Two-dimensional Internal waves in a continu- 
ously stratified fluid 

[References]: 

C.S. Yih, 1965, Dynamics of Inhomogeneous Fluids, MacMillan. 

O. M. Phillips, 1977, Dynamics of the Upper Ocean, Cambridge U. Press. 

P. G. Baines, 1995, Topographical Effects in Stratified Flows Cambridge U. Press. 

M. J. Lighthill 1978, Waves in Fluids , Cambridge University Press. 

Due to seasonal changes of temperature, the density of water or atmosphere can have 
significant variations in the vertical direction. Variation of salt content can also lead to 
density stratification. Freshwater from rivers can rest on top of the sea water. Due to 
the small diffusivity, the density contrast remains for a long time. 

Consider a calm and stratified fluid with a static density distribution p (z) which 
decreases with height (z). If a fluid parcel is moved from the level z upward to z + (, it 
is surrounded by lighter fluid of density p(z + dz). The upward buoyancy force per unit 
volume is 

d~o 

g(p(z + C) - p{z)) » g—C 

dz 

and is negative. Applying Newton's law to the fluid parcel of unit volume 

_d 2 ( dp 
^ dt 2 ^ dz 

or 

§ + N°C = (8.1) 

where 

is called the Brunt-Vaisala frequency. This elementary consideration shows that once 
a fluid is displaced from its equilibrium position, gravity and density gradient provides 
restoring force to enable oscillations. In general there must be horizontal nonunifomities, 
hence waves are possible. 

We start from the exact equations for an inviscid and incompressible fluid with 
variable density. 
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For an incompressible fluid the density remains constant as the fluid moves, 

pt + q • Vp = (8.3) 

where q = (u,w) is the velocity vector in the vertical plane of (x,z). Conservation of 
mass requires that 

V-q = (8.4) 
The law of momentum conservation reads 

p(qt + q • Vq) = - Vp - pge z (8.5) 
and e z is the unit vector in the upward vertical direction. 

8.1 Linearized equations 

Consider small disturbances 

p = p + p', p = p{z) + p', q=(u',w') (8.6) 

with 

p > p', p > p' (8.7) 

and u', v', w' are small. Linearizing by omitting quadratically small terms associated 
with the fluid motion, we get 

P't + = 0. (8.8) 

u' x + w' z = (8.9) 

pu't = -P' x (8-10) 

pw't = ~Pz ~ P'z -9P- 9p' (8.H) 
In the last equation the static part must be in balance 

= ~P Z - 9P, (8-12) 

hence 

p(z) = / p(z)dz. (8.13) 
Jo 
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The remaining dynamically part must satisfy 



Piv' t = -P'z-9P' (8-14) 



Upon eliminating p 1 from the two momentum equations we get 



dp 

—u' t + p(u' z - w' x ) t = gp' x (8.15) 
Eliminating p' from (8.8) and (8.15) we get 

j-v! tt + p(u' z - w' x ) tt = gp' xt = -9j- z w 'x (8-16) 
Let us introduce the disturbance stream function tp: 

v! = fa, w' = -fa (8.17) 

It follows from (8.16) that 

d — 

P (4>xx + tpzz) tt = (94>xx - Ipztt) (8-18) 



by virture of Eqns. (8.8) and (8.17). Note that 



is the Brunt-Vaisala frequency. In the ocean, density gradient is usually very small ( 
iV ~ 5 x 10~ 3 rad/sec). Hence p can be approximated by a constant reference value, say, 
Po = p(0) in (8.10) and (8.14) without much error in the inertia terms. However density 
variation must be kept in the buoyancy term associated with gravity, which is the only 
restoring force responsible for wave motion. This is called the Boussinesq approximation 
and amounts to taking p to be constant in (Eq:17.1) only. With it (8.18) reduces to 

^xx + ^zz) tt + N 2 (z)ij xx = 0. (8.20) 

Note that because of linearity, u' and w' satisfy Eqn. (8.20) also, i.e., 

Kx + Kz)tt + N 2 w' xx = (8.21) 

etc. 
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8.2 Linearized Boundary conditions on the sea surface 

Dynamic boundary condition : Total pressure is equal to the atmospheric pressure 

(p + p') z=( = 0. (8.22) 

On the free surface z = C, we have 



P~-9 [ p(0)dz = -gp(0)C 
Jo 



Therefore, 



implying 



Kinematic condition 



-pg( + p' = 0, z = 0, (8.23) 

-pgQxt = -p' xxt , z = 0. (8.24) 



Ct = w, z = 0. (8.25) 
The left-hand-side of (8.24) can be written as 

-pgCxxt = -pg w' xx 

Using 8.10, the right-hand-side of 8.24 may be written, 

-Pxxt = pU x tt = -P w 'ztt 

hence 

w 'ztt - 9 w 'xx = °> 011 2 = 0- (8.26) 
Since w' = —ip x , ip also satisfies the same boundary condition 

ipztt - gipxx = 0, on z = 0. (8.27) 

On the seabed, z = —h(x) the normal velocity vanishes. For a horizontal bottom we 
have 

ip(x,-h,t) = 0. (8.28) 
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Figure by MIT OCW. 

Figure 20: Typical variation of Brunt-Vaisala frequency in the ocean. From O. M. 
Phillips, 1977 

9 Internal waves modes for finite N 



Consider a horizontally propagating wave beneath the sea surface. Let 



±ikx „— iujt 



tf> = F(z) e ±lkx e 



(9.1) 



From Eqn. (8.2F 



-co 



or, 



dz 2 

N 2 -uj 2 

+ 5— 



dz 2 



k 2 F = z < 0. 



On the (horizontal) sea bottom 



F = z = -h. 



(9.2) 



(9.3) 



From Eqn. (8.27), 



— -g — F = z = 0. 



(9.4) 



dF e 

dz ^ LO 1 

Equations (9.2), (9.3) and (9.4) constitute an eigenvalue condition. 

If lo 2 < N 2 , then F is oscillatory in z within the thermocline. Away from the 
thermocline, lo 2 > N 2 , W must decay exponentially. Therefore, the thermocline is a 
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waveguide within which waves are trapped. Waves that have the greatest amplitude 
beneath the free surface is called internal waves. 

Since for internal waves, uj < N while iV is very small in oceans, oceanic internal 
waves have very low natural frequencies. For most wavelengths of practical interests 
uj 2 <C gk so that 

F = on z = 0. (9.5) 

This is called the rigid lid approximation, which will be adopted in the following. 

With the rigid-lid approximation, if iV=constant (if the total depth is relataively 
small compared to the vertical scale of stratification, the solution for F is 



\/N 2 - UJ 2 

F = Asm [ k(z + h) - ) (9.6) 



where 



JN 2 - uj 2 

kh- = nir, n = 1,2,3... (9.7) 

UJ 

This is an eigen-value condition. For a fixed wave number k, it gives the eigen-frequencies, 

N , s 

uj n = (9.8) 

V 1 + (ff ) 

For a given wavenumber k, this dispersion relation gives the eigen-frequency co n . For a 
given frequency uj, it gives the eigen-wavenumbers k n , 

k n = ~, ; (9.9) 

h ViV 2 - uj 2 V ; 

For a simple lake with vertical banks and length L, < x < L, we must impose the 
conditions : 

u' = 0, hence ip = 0, x = 0, L (9.10) 

The solution is 



tp = Asm k m x exp(— iuj nm t) sin 

with 



k m (z + h) 



UJ 



2 



UJ-r 



(9.11) 



k m L = rmr, m = l,2, 3, ... (9-12) 

The eigen-frequencies are: 

N 

w m = , = (9.13) 

V 1 + (S) 2 
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10 Internal waves in a vertically unbounded fluid 

Consider N = constant (which is good if attention is limited to a small vertical extent), 
and denote by (a, (3) the (x, z) components of the wave number vector k Let the solution 
be a plane wave in the vertical plane 

Then 



to 2 = N 2 



2 



tt 2 + /? 2 

or 



lo = ±N^ (10.1) 

k 

k 2 = a 2 + ft 2 (10.2) 
This is the dispersion relation. Note that is 

-^ = ±cos#' (10.3) 

where 9' is the inclination of k with repect to the x axis. For a given frequency, there 

are two possible signs for a. Since the above relation is also even in j3, there are four 
possible inclinations for the wave crests and troughs with respect to the horizon; the 
angle of inclination is 

\9'\ = cos" 1 -^ (10.4) 

For d > N, there is no wave. 

To under the physics better we note first that the phase velocity is 

— * 

d = U = ± ¥ M) (10 - 5) 



while the group velocity components are 



(10.6) 





a 






a- 










¥ 


k) 





ygx da J ~" \k k 2 k ) k V k 2 J ^k 3 ' 



Cgz ~ dp " T F- (1 °- 7) 



Thus 
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Cg 




Figure 21: Phase and group velocities 

Two conclusions can be drawn. First, the group velocity is perpendicular to the phase 
velocity, 

C 9 -C = 0. (10.9) 

Second, 

C + C g = ±^ (a 2 + (3 2 , 0) = ±^ (k, 0) (10.10) 

The vector sum of C and C g is a horizontal vector, as shown by any of the sketches in 
Figure 17. Note that when the phase velocity as an upward component, the group veloc- 
ity has a downward component, and vice versa. Now let us consider energy transport, 
from (8.10) we get 

-p' x = pt{; zt = puj^ e^ + ^-^ 
hence the dynamic pressure is 

p' = iuSp-iioe 1 ^^-^ (10.11) 
a 

The fluid velocity is easily calculated 

4 = („', „') = (^, -^) = iptf, -a)^ ax+ ^-^ (10.12) 
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The averaged rate of energy transport is therefore 

E = \f\^\ 2 ^{^-a) (10.13) 

which is in the same direction of the group velocity. 

Energy must radiate outward from the oscillating source, hence the group velocity 
vectors must all be outward. Since there are 4 directions for k. There are four radial 
beams parallel to S g , in four quadrants, forming St. Andrews Cross. The crests (phase 
lines) in the beam in the first quadrant must be in the south-easterly direction. Similarly 
the crests in all four beams must be outward and toward the horizontal axis. Let 9 be 
the inclination of a beam (i.e., (C g ) with respect to the y axis, then 9 = tt/2 — 9' in the 
first quadrant. The dispersion relation can be written as 

-^ = ±sin# (10.14) 

where 9 is the inclination of a beam and not of the wavenumber vector. 

Movie records indeed confirm these predictions. Within each of the four beams which 
have widths comparable to the cylinder diameter, only one or two wave lengths can be 
seen. 

This unique property of anisotropy has been verified in dramatic experiments by 
Mowbray and Stevenson. By oscillating a long cylinder at various frequencies verti- 
cally in a stratified fluid, equal phase lines are only found along four beams forming St 
Andrew's Cross, see Figure (22) for uj/N = 0.7, 0.9 and 1.11. It can be verified that 
angles are \9\ = 45° for uj/N = 0.7, and \9\ = 64° for u/N = 0.9. In the last photo, 
lo/N = 1.11. There is no wave. These results are all in accord with the condition (10.4). 

Comparison between measured and predicted angles is plotted in Figure (23) for a 
wide range of uj/N 



11 Reflection of internal waves at a plane boundary 

For another interesting feature, consider the reflection of an internal wave from a slope. 

Recall that 9' = ±cos _1 i.e., for a fixed frequency there are only two allowable 
directions with respect to the horizon. Relative to the sloping bottom inclined at 9 the 
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Figure 22: St Andrew's Cross in a stratified fluid. Top: uj/N = 0.7; bottom left 
u)/N = 0.9; bottom right: lo/N = 1.11. From Mowbray & Rarity, 1965, JFM 
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(0/N 



Figure by MIT OCW. 

Figure 23: Comparison of measured and predicted angles of internal-wave beams. lo/N 
vs. sine?. From Mowbray & Rarity, 1965, JFM 
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inclinations of the incident and reflected waves must be different, and are respectively 
6' + 9 and 9' - 9 , see Figure 24. 

Let £ be along, and i] be normal to the slope. Since the slope must be a streamline, 
ipi + ip r must vanish along r\ = and be proportional to e*^ - ^; the total stream 
function must be of the form 

In particular the wavenumber component along the slope must be equal, 

kf = k { ; ] = a 

Consider first 9 > 9', as shown in the left diagram in Figure(24). We must have 

k (i) cos(9' + 9 ) = k (r) cos(0' - 9 ), 

which implies first that 

k^^k^ r \ (11.1) 

The incident wave and the reflected wave have different wavelengths! Furthermore, the 
reflected wave is directed up-slope.If, however, 9 > 9', as shown in the right diagram in 
Figure(24), then 

JfcW cos(6»' + 9 ) = k^ cos(^ - (9 - 9')) = -JfeW cos(6» - 9'), 

The reflected wave is directed down-slope instead (Phillips, 1977). 
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Figure 24: Internal wave reflected by in inclined surface. Top: 9 < 9'. Bottom: 9 > 9'. 
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